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A NEW CLASSIFICATION METHOD FOR GAMMA-RAY BURSTS 
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ABSTRACT 

Recent 5wra/f observations suggest that the traditional long vs. short GRB classification scheme does 
not always associate GRBs to the two physically motivated model types, i.e. Type II (massive star 
origin) vs. Type I (compact star origin). We propose a new phenomenological classification method 

5 /3 

of GRBs by introducing a new parameter e — E~f^ so ,52/ E p z 2 , where -E 7) i S o is the isotropic gamma-ray 
energy (in units of I0 52 erg), and E PjZ is the cosmic rest frame spectral peak energy (in units of 100 
keV). For those short GRBs with "extended emission", both quantities are defined for the short/hard 
spike only. With the current complete sample of GRBs with redshift and E p measurements, the 
e parameter shows a clear bimodal distribution with a separation at e ~ 0.03. The high-e region 
encloses the typical long GRBs with high-luminosity, some high-z "rest-frame-short" GRBs (such as 
GRB 090423 and GRB 080913), as well as some high-z short GRBs (such as GRB 090426). All 
these GRBs have been claimed to be of the Type II origin based on other observational properties 
in the literature. All the GRBs that are argued to be of the Type I origin are found to be clustered 
in the low-e region. They can be separated from some nearby low-luminosity long GRBs (in 3c) 
by an additional Tgo criterion, i.e. Tqq iZ < 5 s in the Swift/BAT band. We suggest that this new 
classification scheme can better match the physically-motivated Type II/I classification scheme. 
Subject headings: gamma-ray bursts: general — methods: statistical 



1. INTRODUCTION 

Phenomenologically, gamma-ray bursts (GRBs) are 
classified as long vs. short with a division line at the 
observed duration T 90 ~ 2 s (Kouveliotou ct al. 1993). 
Robust associations of the underlying supernovae (SNe) 
with some long GRBs (Galama et al. 1998; Stanek et al. 
2003; Hjorth et al. 2003; Malesani et al. 2004; Modjaz 
et al. 2006; Pian et al. 2006) and the fact that long 
GRB host galaxies are typically irregular galaxies with 
intense star formation (Fruchter et al. 2006) suggest that 
they are likely related to the deaths of massive stars, 
and the "collapsar" model has been widely recognized 
as the standard scenario for long GRBs (Woosley 1993; 
Paczyhski 1998; Woosley & Bloom 2006). Observational 
breakthroughs led by the Swift mission (Gehrels et al. 
2004) suggest that at least some short GRBs are associ- 
ated with nearby host galaxies with little star formation 
(Tanvir et al. 2005), and that they are not associated 
with an underlying SN (Gehrels et al. 2005; Villasenor 
et al. 2005; Fox et al. 2005; Hjorth et al. 2005; Berger 
et al. 2005), favoring the idea that they are produced by 
mergers of two compact stellar objects, such as NS-NS 
and NS-BH mergers (Eichler et al. 1989; Narayan et al. 
1992; Nakar 2007). 

However, several lines of observational evidence in the 
Swift era suggest that duration is not necessarily a re- 
liable indicator of the physical nature of a GRB. (1) 
The non-detection of a SN signature associated with the 
nearby long GRBs 060614 and 060505 (Gehrels et al. 
2006; Gal- Yam et al. 2006; Fynbo et al. 2006; Delia 
Valle et al. 2006) disfavors the conventional collapsar 
scenario of long GRBs (Woosley & Bloom 2006). Some 
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observational properties of the ~ 100 s long GRB 060614 
are very similar to those of some short GRBs (Gehrels et 
al. 2006; Zhang et al. 2007), making it likely associated 
with the physical category that most short GRBs belong 
to. (2) Significant soft "extended" gamma-ray emission 
and late X-ray flares are observed in a handful of "short" 
GRBs (Barthelmy et al. 2005; Norris & Bonnell 2006; 
Per ley et al. 2009), suggesting that they are not neces- 
sarily short. A new physically motivated classification 
scheme, i.e., Type II (massive star origin) vs. Type I 
(compact star origin) was proposed (Zhang 2006; Zhang 
et al. 2007). Kann et al (2010, 2008) systematically stud- 
ied the optical afterglow emission properties of Type II 
and Type I GRBs, and suggested that these properties 
carry important information about the nature of GRB 
progenitors. A full definition of Type I /II GRBs as well 
as the multiple observational criteria are presented in 
Zhang et al. (2009). (3) Two high-redshift gamma-ray 
bursts, GRB 080913 at z = 6.7 (Greiner et al. 2009; 
Perez-Ramirez et al. 2010) and GRB 090423 at z = 8.2 
(Tanvir et al. 2009; Salvaterra et al. 2009), appear as 
intrinsically short GRBs. Their observed durations are 
(8±l)s and ~ 10.3s in the Swift/BAT band, respectively. 
Corrected redshift, the durations of two high-z GRBs arc 
shorter than 2 seconds in the rest frame. The physical 
origin of these high-z GRBs have been subject to debate 
(Greiner et al. 2009; Perez-Ramirez er al. 2010; Tanvir 
et al. 2009; Salvaterra et al.2009). Although compact 
star mergers may occur at such a high redshift (Belczyn- 
ski et al. 2010), various observational properties of these 
two GRBs point towards the massive star origin (Zhang 
et al. 2009; Lin et al. 2009; Belczynki et al. 2010). 
(4) A more striking case is GRB 090426, whose observed 
BAT band Tgo is only 1.2 ± 0.3s, and the rest frame du- 
ration is only ~ 0.33 s at z = 2.609 (Levesque et al. 
2009). It greatly exceeds the previous short GRB red- 
shift record, i.e. z = 0.923 for GRB 070714B (Graham et 
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al. 2009). Phenomenologically, this is an unambiguous 
short-duration GRB, but the available afterglow and host 
galaxy properties point towards a different picture: it has 
a blue, very luminous, star-forming putative host galaxy 
with a small angular offset of the afterglow location from 
the center, and a similar medium density as typical Type 
II GRBs. All these suggest that the burst is more closely 
related to the core collapse of a massive star (Levesque 
et al. 2009; Antonelli et al. 2009; Xin et al. 2010). (5) 
Monte Carlo simulations suggest that the compact star 
merger model cannot interpret both the Swift z-known 
short GRB sample and the BATSE short GRB sample. 
Instead one may need a mix of GRBs from compact star 
mergers and from massive star core collapses to account 
for the observed short GRB population (Virgili et al. 
2009; Cui et al. 2010). 

In summary, the traditional long vs. short GRB 
scheme no longer always correspondingly match the two 
distinct physical origins, i.e., collapse of massive stars 
(Type II GRBs) vs. mergers of compact stars (Type I 
GRBs) (Zhang 2006; Zhang et al. 2007, 2009; Bloom et 
al. 2008). Zhang et al (2009) proposed a procedure (a 
flow chart, see their Fig. 8) invoking a full set of observa- 
tional characteristics, including supernova(SN)-GRB as- 
sociation, specific star-forming rate (SFR) of the host 
galaxy, burst location offset, burst duration, hardness, 
spectral lag, statistical correlations, energetics, and af- 
terglow properties, to judge the physical category of a 
GRB. On the other hand, for most GRBs one may not 
immediately retrieve all the information needed for such 
a classification. It would be interesting to search for ad- 
ditional phenomenological classification schemes to see 
whether there exist other quantities, especially those in- 
voking prompt GRB emission properties only, that can 
give a good indication of the physical nature of a GRB. 
This motivates us to explore a new phenomenological 
classification method. Besides the burst duration, we 
think both the burst energy and the spectral properties 
in the burst rest frame would carry important informa- 
tion about the nature of the burst. In this paper, we 
propose a new discriminator based on the isotropic burst 
energy and rest frame peak energy (E p ) of the vf u spec- 
trum of the prompt gamma-rays. We will show that this 
classification method is more closely connected to the 
Type II /I physical classification scheme. 

2. SAMPLE 

In order to develop our classification method, we use a 
sample of all GRBs with both redshift and spectral pa- 
rameters known up to Mar, 2010. Our sample includes 
137 GRBs. They are detected by BeppoSAX, HETE- 
2, Swift, Suzaku, and Fermi. Most of them are taken 
from Amati et al. (2008, 2009), Kann et al. (2010), 
Krimm et al. (2009), Zhang et al. (2009) and references 
therein. Their isotropic gamma-ray energy (E ly i SO ) is 
calculated in the GRB rest-frame 1 — 10 4 kev band with 
the spectral parameters in order to avoid instrument se- 
lection effect. The soft XRF 080109 (Soderberg et al. 
2008), whose emission is in the XRT band (0.3-10kev) 
instead of BAT band, is also included in our sample. 
The XRT and UVOT observations place a limit of its 
E p in the range 0.037 kev < E pea k < 0.3 kev. We 

adopt E peak » 0.12±2:os9 kev ( Li 2008 )- Ei S ht bursts 



in our sample (GRBs 050709, 050724, 051210, 060614, 
061006, 061210, 070714B, and 071227) are the so-called 
short GRBs with "extended emission" discussed in lit- 
erature. For these GRBs, both parameters (E 7 ^ so and 
E p ,z) are derived for the initial hard spike only (extended 
emission excluded)0. We tabulate only those GRBs with 
T 90 , z = T 90 /(l + z) < 2 s (Twenty-nine GRBs) in our 
sample in Table 1. 

3. A NEW GRB CLASSIFICATION PARAMETER 

Firstly, we show the distribution of Tgo. 2 in Fig. 1(a). 
Similar to that shown in Lin et al. (2009), Tgo yZ for the 
current GRB sample with redshift measurements has a 
tentative bimodal distribution, with peaks at logTgo iZ = 
—0.5 and 1.25, respectively. Due to the overlap between 
the two classes (long vs. short), those GRBs having 
2go,z = 1 ~ 5 seconds cannot be unambiguously cate- 
gorized into either group0. We test the bimodality with 
the KMM algorithm presented by Ashman et al. (1994), 
and obtain a chance probability p ~ 10 ~ 2 . The correct 
allocation rate for short GRBs estimated with the KMM 
algorithm is 73.3%. 

We develop our classification method by defining a pa- 
rameter, 

e{k) = £ 7 , iso , 52 /i% Zi2) (1) 

where E lfis o,m = £ 7 ,i so /10 52 ergs and £ P)Z = E p (l + 
z)/10 2 keV. Note that e does not depend on any spec- 
trum - energy correlation or theoretical assumption. We 
explore various n values, and found that given a same 
observed GRB (known fluence and E p ) the implicit z- 
dependence of e(n) cx D 2 L {z)/(l + z) 1+K essentially van- 
ishes at z > 2, if k = 5/3 (see Fig. 2 for a comparison 
of the z-dependence of e for different k). We therefore 
define 

e = E lAso , 52 /E pz2 , (2) 

which removes the redshift dependence for high- z GRBs. 
This parameter may be more suitable for the classifica- 
tion of high-z GRBs, such as GRBs 090423 and 080913. 

We calculate e for the bursts in our sample and show 
the loge distribution in Fig.l (b). A clear bimodal fea- 
ture is found, with a division line at ~ 0.03 (high-e vs. 
low-e). The high-e portion follows a log normal distri- 
bution with e = 0.80 ± 0.11 (lcr). For the low-e it is 
0.003 ± 0.002. We test the bimodality of the e distri- 
bution with the KMM algorithm and obtain a chance 

3 We note that Tqq depends on the detector's sensitivity and 
energy band. In our analysis, we define Tgo using the BAT energy 
band and sensitivity. In the Swift era, it is found that some "short" 
GRBs (as observed in the BATSE band) have extended emission. 
It is the convention of the Swift team and the community to define 
a "short" GRB based on the duration of the short spike only, with 
the extended emission excluded. Observationally, the spectrum of 
the short spike is much harder than that of the extended emission, 
the latter shares many properties with X-ray flares. For a same 
burst, the "extended emission" depends on redshift, since it may 
be buried within the background if the redshift is high enough. 
We therefore only use the information of the short/hard spike to 
perform our classification. 

4 It was proposed that those GRBs having Tqq iZ = 1 ~ 5 sec- 
onds may be an intermediate population based on the excess of 
the GRBs over the bimodal distribution of Tgo, z . Although these 
GRBs are less energetic and have dimmer afterglows than typical 
long GRBs, they share many similar properties of typical GRBs, 
such as spectral lag and spectrum-energy correlations. They may 
be a sub-group of the long GRBs (de Ugarte Postigo 2010). 
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probability p < 1CP 4 . The overall correct allocation rate 
of GRBs in the two categories (high-e vs. low-e) is as 
high as 99.8%, indicating that £ is a new parameter for 
GRB classification. 

It is generally believed that GRBs originate either from 
core collapses of massive stars (Type II) or mergers of 
compact stars (Type I). One may ask whether this phe- 
nomenal classification matches the physical scheme of 
the Type I and Type II. The origins of most GRBs in 
our sample have been extensively discussed in literature 
(Zhang et al. 2009; Kann et al. 2010, 2008). We there- 
fore examine how the members in our two new categories 
(high-e vs. low-e) are associated with the two physi- 
cal model types. To define the physical category of the 
GRBs, we apply the flow chart of Zhang et al. (2009, 
their Fig. 8). We denote all "Type II" or "Type II can- 
didates" as "Type II", and "Type I" or "Type I can- 
didates" as "Type I" (which include the Type I Gold 
sample and most other short/hard GRBs in Zhang et al. 
2009). We note that these definitions are fully consistent 
with those adopted in Kann et al. (2010, 2008). Fig. 1 
displays all the GRBs in our sample in the e — Tgo.z plane. 
Different symbols are used to mark different groups. For 
example, triangles denote the GRBs that are believed to 
be of the Type II origin, among which the black triangles 
denote the traditional long GRBs, and the blue triangles 
denote the GRBs those with Tgo jZ 's shorter than 2 s. The 
red circles denote those GRBs that are believed to have 
the Type I origin. Six green stars denote nearby low 
luminosity GRBs with supernova associations. It is re- 
markable to see that the e classification clearly separates 
Type II GRBs from Type I GRBs. 

In order to quantitatively assess the clustering fea- 
ture among different types of GRBs, we derive the GRB 
distribution probability p(loge, logTgo jZ ) from the ra- 
tio between the number of GRBs in the grid loge + 
dloge, logT 90)Z + G?logTgo, 2 and the total number of 
GRBs in our sample. Since the distribution of e is con- 
cerned in our analysis, we take the error of loge into 
account in our calculation of p. For a given GRB with 
loge ± S loge, we assume that the probability distribu- 
tion of this GRB is a normalized log-normal distribution 
centering at loge with a width of S loge. The proba- 
bility (pi) of the GRB in a given rage [loge — dloge/2, 
loge + dloge/2] is calculated by integrating the proba- 
bility distribution function in the range. The probability 
p(loge, logTgo i2 ) is then obtained by summing up pi over 
the GRBs in our sample. We take c?logTgo !2 = 0.35 
and dloge = 0.60 in our calculation. The contours 
of the probability distribution are shown in Fig. 1(c). 
One can clearly identify two clustering regions [with 
p(\og e, log Tgo jZ ) > 0.1] centered in the high-e regime and 
the low-e, low-Tgo iZ regime. At the 3cr level, the high-e 
group includes typical long GRBs with high-luminosity, 
some intrinsic short duration GRBs (including the high-z 
GRB 090423 and GRB 080913), as well as some high- 
z short GRBs (e.g. GRB 090426). In contrast, the 
low-e group contains the extensively discussed Type I 
GRBs, including those with and without extended emis- 
sion. Based on the probability contours, a low-e GRB 
with T 90 , 2 < 5 s in the Swift/BAT band would be 
identified as a Type I GRB at the 3<r confidence level. 
GRBs 060614 and 060505 are marginally included in the 
p(loge, logTgo, z ) > 0.003 region (the light grey region in 



Fig. lc) of the Type I GRBs. Nearby low luminosity long 
GRBs (LL-GRBs), e.g. GRBs 980425, 031203, 050826, 
and 060218 (green stars in Fig. lc) are out of the 3 a 
contours of the Type I and Type II GRBs. This may hint 
a distinct GRB population as suggested by Liang et al. 
(2007; see also Soderberg et al. 2004; Cobb et al. 2006). 

4. CONCLUSIONS AND DISCUSSION 

By introducing a new parameter e, we have proposed 
a new phenomenological classification method for GRBs. 
We demonstrated that GRBs with known redshifts are 
cozily classified into two classes with a separation e = 
0.03. Due to the clear bimodal distribution of e with lit- 
tle overlap, the overall correct allocation rate of a GRB 
to a particular category (high-e vs. low-e) is as high 
as > 99.8%, which is much better than the traditional 
duration classification method. More importantly, this 
classification scheme is more closely related to the two 
physically motivated model classes - Type II vs. I. In 
particular, the high-e category is a good representation 
of the Type II GRBs, while the low-e category, with an 
additional duration criterion Tgo^ < 5 s, is a good rep- 
resentation of the Type I GRBs. We suggest that the 
e parameter can be evaluated for future GRBs with z 
measurements, and the e-based classification may be per- 
formed regularly to infer the physical origin of a GRB. 

The well-separated bimodal distribution of loge de- 
pends on our selection of k in Eq. (1). As shown in Fig. 
2, e is insensitive to z at z > 2 for a given burst by taking 
k = 5/3, but it is not at low redshift. One may suspect 
that if the bimodal distribution is due to the redshift ef- 
fect. To examine this issue, we re-plot the Figure 1 by 
separating GRBs into the z > 2 and the z < 2 samples 
(Figure 3, left and middle panels, respectively). We also 
show the GRB distributions in the log(l + z) — log(e) 
plane for all the bursts in our sample in Figure 3. It 
is found that the the clustering feature of loge is not 
caused by the redshift selection effect. Therefore, e is a 
GRB discriminator for both GRBs at z < 2 and z > 2. 

Another concern about our classification method may 
be sample uniformity. The GRBs in our sample were 
detected with different instruments with different instru- 
mental threshold. However, our classification method 
is based on the global energy and spectral information. 
Different from the burst duration, a well-measured E p 
essentially does not depend on the detector. So is E^^ so . 
The E p values of the GRBs in our sample, which expand 
almost four orders of magnitude, are well within the en- 
ergy band of previous and current GRB missions, so they 
are well measured. The E- lso values expand almost eight 
orders of magnitude. Therefore, our sample could be re- 
garded as a complete sample with known E p and E mo for 
current missions. 

We should note that our classification scheme is par- 
ticularly helpful to diagnose the physical origin of GRBs 
with a rest-frame short duration (e.g. Tgo j2 < 2 s, Table 
1). For example, it convincingly classifies GRBs 090423, 
080913 and 090426 into the high-e, and therefore, the 
Type II category, which is otherwise difficult with the du- 
ration criterion^]. The distinction of the Tqq,z < 2s bursts 

5 As shown in Figure 2, a GRB with the same observed fluence 
and E p at higher redshift may have a higher e. However, the 
bimodal distribution may not be due to the redshift effect (see 
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with different e values is also evident from their X-ray 
afterglow properties. Kann et al. (2008) compared the 
optical afterglows Type I GRBs with Type II GRBs, and 
found that those of Type I GRBs have a lower average 
luminosity and show a larger intrinsic spread of luminosi- 
ties. In Figure 4, we compare the XRT lightcurves and 
the 12-hour X-ray luminosities between the two types of 
T 90 , z < 5s s GRBs. The low-e T go , z < 2s s GRBs (Type 
I, red) are systematically fainter than the high-e ones 
(Type II, bluefl whose properties are rather similar to 
other high-e long GRBs (black). 

Physically it is not obvious why the parameter e (cor- 
responding to k = 5/3) gives a cozy classification scheme. 
We have tried other k values in Eq.(l), but the classifi- 
cations defined by other e(n) are not as clean as the one 
defined by Eq.(2). As shown in Figure 2, an apparent 
advantage of k = 5/3 is to diminish the z-dependence at 
high- z for a given observed GRB. 

One can also discuss the connection between the log- 
normal distribution of e in the high-e regime and some 
empirical relations of Type II GRBs. With the Amati- 
relation E p , z oc E°-f* (Amati et al. 2009), wc find that 

e oc -Ep^ 85 , insensitive to _E PiZ . Liang & Zhang (2005) 
discovered an empirical relation among -E 7 .iso, E PiZ , and 
t b , z , namely, £ 7 , iso oc ±oa \*- 24±oi3 , where t b , z 
is the rest-frame break time of the optical afterglow 
lightcurve. This gives e oc J E°; z 7±017 < b " z - 24±0 - 23 , which 
is sensitive to tb jZ . A clustering in e therefore corre- 
sponds to a clustering in tb,z- Interpreting the £t>,z as the 
jet break time, the Liang-Zhang relation may be trans- 



lated into the Ghirlanda-relation, £' 7! iso.52(l — cosOj) ~ 
AE~IJ, where A ~ 0.01 - 0.03 (Ghirlanda et al. 2004; 
Dai et al. 2004). Given e = E^^o^/ E p l 2 , we can get 

e w A(l — costfj) -1 , where 8j is the jet opening angle. 
The lognormal e distribution (e = 0.80 ± 0.11, at ler) 
therefore corresponds to a lognormal jet opening angle 
distribution for Type II GRBs, i.e., 9j ~ 0.16 ± 0.03 rad 
(assuming A ~ 0.01). We note that the jet opening an- 
gle derived by Ghirlanda et al. (2004) indeed clusters in 
a small range. Since the Ghirlanda-relation is very dif- 
ficult to understand physically, we suspect that it may 
be related the fact that e is quasi-universal for Type II 
GRBs. Future data with very early or very late opti- 
cal break times may test whether the Ghirlanda/Liang- 
Zhang relations or the fact of a quasi-universal s are more 
fundamental. 

We acknowledge the use of the public data from the 
Swift and Fermi data and GCN circulars archive. We 
thank helpful comments from the referee. This work is 
supported by the National Natural Science Foundation 
of China (Grants 11025313, 10873002, 11078008), the 
National Basic Research Program ("973" Program) of 
China (Grant 2009CB824800), Guangxi SHI-BAI-QIAN 
project (Grant 2007201), Guangxi Science Foundation 
(2010GXNSFC013011), the program for 100 Young and 
Middle-aged Disciplinary Leaders in Guangxi Higher 
Education Institutions, and the research foundation of 
Guangxi University. It is also partially supported by 
NASA NNX09AT66G, NNX10AD48G, and NSF AST- 
0908362. 



REFERENCES 



Abdo, A. A., et al. 2009, Nature, 462, 331 

Amati, L., Guidorzi, C, Frontera, F., et al. 2008, MNRAS, 391, 
577 

Amati, L., Frontera, F., & Guidorzi, C. 2009, A&A, 508, 173 
Antonelli, L. A., et al. 2009, A&A, 507, L45 
Ashman, K. M., Bird, C. M., & Zepf, S. E. 1994, AJ, 108, 2348 
Barthelmy, S. D., et al. 2005, Nature, 438, 994 
Belczynski, K., Holz, D. E., Fryer, C. L., et al. 2010, ApJ, 708, 117 
Berger, E., Price, PA., Cenko, S.B., et al. 2005, Nature, 438, 988 
Bloom, J. S., Butler, N. R., & Perley, D. A. 2008, AIPC, 1000, 11 
Butler, N. R., Kocevski, D., Bloom, J. S., et al. 2007, ApJ, 671, 
656 

Cenko, S. B., et al. 2010, GCN, 10389, 1 

Cobb, B. E., Bailyn, C. D., van Dokkum, P. G, &; Natarajan, P. 

2006, ApJ, 645, L113 
C ui, X.-H., Aoi, J Nagataki, S., & Xu, R.-X. 2010, 

larXiv:1004.2302l 
Dai, Z. G., Liang, E. W., & Xu, D. 2004, ApJ, 612 L101 
de Ugarte Postigo, A., et al. 2010. larXiv: 1006.44691 
Delia Valle, M., et al. 2006, Nature, 444, 1050 

Eichler, D., Livio, M., Piran, T., & Schramm, D. N. 1989, Nature, 
340, 126 

Fox, D.B., Frail, DA., Price, P.A., et al. 2005, Nature, 437, 845 
Fruchter, A.S., Levan, A. J., Strolger, L., et al. 2006, Nature, 441, 
463 

Fynbo, J.P.U., Watson, D., Thone, C.C., et al. 2006, Nature, 444, 
1047 

Galama, T. J., et al. 1998, Nature, 395, 670 

Gal- Yam, A., Fox, D.B., Price, PA., et al. 2006, Nature, 444, 1053 
Gehrels, N., et al. 2004, ApJ, 611, 1005 

Figure 3). 

6 We cannot exclude the possibility that this is a selection effect, 
since high-z Type I GRBs may have fainter X-ray afterglows and 
therefore may not be detected. 



Gehrels, N., Sarazin, C.L., O'Brien, P.T., et al. 2005, Nature, 437, 
851 

Gehrels, N., Norris, J. P., Barthelmy, S.D., et al. 2006, Nature, 
444, 1044 

Ghirlanda, G, Ghisellini, G, & Lazzati, D. 2004, ApJ, 616, 331 
Graham, J. F., et al. 2009, ApJ, 698, 1620 

Greiner, J., Kriihler, T., McBreen, S., et al. 2009, ApJ, 693, 1610 
Hjorth, J., Sollerman, J., M0ller, P., et al. 2003, Nature, 423, 847 
Hjorth, J., et al. 2005, N ature, 437, 859 
Kann, D. A., et al. 2008, arXiv:0804.1959 
Kann, D. A., et al. 2010, ApJ, 720, 1513 

Kouveliotou, C, Meegan, C. A., Fishman, G. J., et al 1993, ApJ, 
413, L101 

Krimm, H. A., et al. 2009, ApJ, 704, 1405 
Li, L.-X. 2008, MNRAS, 388, 603 
Levesque, E. M., et al. 2009, MNRAS, 1657 
Liang, E.-W., & Zhang, B. 2005, ApJ, 633, 611 
Liang, E. -W, Zhang, B., Virgili, F., & Dai, Z. G. 2007, ApJ, 662, 
1111 

Liang, E.-W., Lii, H.-J., Hou, S.-J., et al. 2009, ApJ, 707, 328 
Lin, L., Liang, E. -W., & Zhang, S. N. 2010, Science in China G: 

Physics and Astronomy, 53, 64 
Malesani, D., et al. 2004, ApJ, 609, L5 
Modjaz, M., et al. 2006, ApJ, 645, L21 
Nakar, E. 2007, Phys. Rep., 442, 166 

Narayan, R., Paczynski, B., & Piran, T. 1992, ApJ, 395, L83 
Norris, J. P., & Bonnell, J. T. 2006, ApJ, 643, 266 
Paczynski, B. 1998, ApJ, 494, L45 
Perez-Ramirez, D., et al. 2010, A&A, 510, A105 
Perley, DA., Metzger, B.D., Granot, J., et al. 2009, ApJ, 696, 
1871 

Pian, E., Mazzali, PA., Masetti, N., et al. 2006, Nature, 442, 1011 
Piranomonte, S., et al. 2006, GCN, 5626, 1 
Sakamoto, T., et al. 2008, GCN, 7761, 1 
Salvaterra, R., et al. 2009, Nature, 461, 1258 



A New Classification Method For GRBs 



5 



TABLE 1 

The sample of 29 GRBs for Tgo^ < 2 s in the rest frame. 



GRB 


z 


Tbo,2 


EE 


E p 


R ■ 
iso 

(10 52 erg) 


log e 


Type 


Afterglow 


Refs. 


name 


redshift 


(s) 


(s) 


(KeV) 






X/0(IR) 




050406 


2.44 


1.57 


n/a 


25+H 


14+0- 13 
u,i — 03 


-0.74+g-™ 


high-e (II) 


Y / Y 


(1) 


050416A 


0.6535 


1.51 


n/a 




o 1+0.01 
u- 1 — 01 


0.00l°;} 2 


high-e (II) 


Y / Y 


(2) 


050509B 


0.2248 


0.06 


n/a 


° Z -80 


2.412' 4 x 10- 4 


-3.62lJ- 2 i 


low-e (I) 


N / N 


(1,3) 


050709 


0.1606 


0.09 


130±7 


83+H 


2.7+" x 10~ 3 


-2.54±g-i5 


low-e (I) 


N / N 


(1.3) 


050724* 


0.2576 


2.39 


154±1 


iio+2S° 


9+_l 9 x 10- 3 


9 9 o + 0.74 
30 


low-e (I) 


Y / N 


(1,3) 


050813 


~0.72 


~0.26 


n/a 


210I™ 


l-5±g | x lO" 2 


9 7=4-1-12 
z - ' °_o 45 


low-e (I) 


N / N 


(1,3) 


050922C 


2.198 


1.41 


n/a 


i3oiijj 




-0.3ll°'il 


high-e (II) 


Y / Y 


(2) 


051221A 


0.5464 


0.91 


n/a 




0.28+°f 


1 00+0.13 
— l-o»_o 17 


low-e (I) 


Y / N 


(3) 


060206 


4.048 


1.51 


n/a 


78±g 


4.318:? 


-0.36l°i° 


high-e (II) 


Y / Y 


(2) 


060502B 


0.287 


0.10 


n/a 


340+™ 


3+2 x 10" 3 


_n CQ+1.14 

o.oy_Q 41 


low-e (I) 


N / N 


(1,3) 


060614* 


0.1254 


~4.44 


106±3 


302±|^ 4 


94+0-04 
u - z ^_ 04 


-1.68+g-g 


low-e (I) 


Y / Y 


(2,3) 


060801 


1.131 


0.23 


n/a 


6201^™ 


n 17 +0.021 
u -! ' _0 021 




low-e (I) 


Y / N 


(1,3) 


060926 


3.208 


1.90 


n/a 


19ti8 


4 2 +0.59 


n i6 +0 - 31 


high-e (II) 


Y / N 


(1.4) 


061006 


0.4377 


~0.35 


120 


640I 227 


0.221°^ 


u -0.26 


low-e (I) 


Y / N 


(1,3) 


061201 


0.111 


0.54 


n/a 


o 7 g+458 
o'J_284 


0.018ig'gg 


o oq+0.38 
, J-'- ,y _0 24 


low-e (I) 


Y / N 


(1,3) 


061210 


0.4095 


0.57 


85±5 


5401™ 


0.09l°; I 


9 r-i+0.82 
*-0±_o 42 


low-e (I) 


N / N 


(1,3) 


061217 


0.827 


0.22 


n/a 


4001"° 


0.031°^ 


-2.96l°' 47 


low-e (I) 


N / N 


(1,3) 


070429B 


0.9023 


0.25 


n/a 


1201™ 


o-o3±g:8i 


9 in+0.75 
z - lo _0 40 


low-e (I) 


N / N 


(1,3) 


070506 


2.31 


1.30 


n/a 


7ll 9 2 5 5 


26+ - 17 


-1 20+ 067 
1-^U_0 26 


high-e (II) 


Y / N 


(1,5) 


070714B 


0.9225 


~1.56 


~100 


1120111° 


1 16+ - 41 
i - iD -0.22 


-2 16+ ' 51 
z - iD -0.25 


low-e (I) 


Y / Y 


(1,3) 


070724A 


0.457 


0.27 


n/a 


~68 


n03+ ° 01 


9 +0.24 
z - J1 -0.24 


low-e (I) 


Y / N 


(3) 


071020 


2.142 


1.34 


n/a 


323+gi 


y '°-4.3 


-0 7n+°' n 
u -' u -0.11 


high-e (II) 


Y / N 


(2) 


071227 


0.394 


1.30 


~100 


100011°° 


n 99 +0.08 
u - zz -0.08 


-2 56+ 014 


low-e (I) 


Y / N 


(3) 


080520 


1.545 


1.10 


n/a 


~30 


n7+°- 02 


-0 94+ 16 
u ' y4 -0.16 


high-e (II) 


Y / N 


(6) 


080913 


6.7 


1.04 


n/a 


1 91 +232 
lzl -39 


7 +1.81 
'—1.81 


_ n 77 +1.39 
u - ' '-0.23 


high-e (II) 


Y / Y 


(4) 


090423 


8.1 


1.13 


n/a 


54+22 
04 -22 


10+3 


-0 11+ - 30 
<J-ll_o.30 


high-e (II) 


Y / Y 


(7,8) 


090426 


2.609 


0.33 


n/a 




42+0- 59 
u -+^-0.04 


n 79 +0.92 
u -' z -0.69 


high-e (II) 


Y / Y 


(9) 


090510 


0.903 


0.26 


n/a 


4414+420 


3 8+ 2 ^ 


-2 fi^+O- 07 


low-e (I) 


Y / N 


(10) 


100206 


0.41 


0.12 


n/a 




6T5+|f x 10~ 2 


-2 61 + 007 
z - Di -0.07 


low-e (I) 


Y / N 


(11) 



al. (2008); (3) Zhang et al. (2009); (4) Piranomonte et al. (2006); (5) Thoene 
(2009); (8) Salvaterra et al. (2009); (9) Levesque et al. (2009); (10) Abdo et 
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Fig. 1. — The 1-D and 2-D distributions for the bursts in our sample in the logTgOjZ vs loge plane along with the bimodal fits (solid 
lines, in panel a, b) as well as the 2-D distribution for the GRBs in our sample in the loge — logTgo, z and loge — log(l + z) plane (panel 
c). The triangles denote the Type II GRB candidates as discussed in literature. Among them the traditional long GRBs are denoted as 
black and the GRBs with Tgo. 2 < 2 s are denoted as blue. The red solid (open) circles represent the Type I GRBs with (without) extended 
emission as discussed in the literature. The diamond denotes the special burst GRB 060505. Nearby LL-GRBs are denoted as green stars. 
The possibility contours for GRB clustering are marked, with p (loge, log Tgo, z ) > 0.1 (dark grey) and p (loge, log Tgo, z ) > 0.003 (light 
grey), respectively. The dashed vertical and horizontal lines are the divisions of e = 0.03 and Tgo, z = 5 seconds, respectively. 



















k= 5/3 , 






fe= 4/3 






k=2 



Fig. 2. — Dependence of e(re) on (1 + z) for k = 4/3, 5/3, 2.0 for a given burst with the same measured fluence and E p . 
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Fig. 3. — The 1-D (panels a and b) and 2-D (panel c) distributions for the bursts in our sample in the logTgo jZ — loge plane for GRBs 
for z > 2 (left group of panels) and z < 2 (middle group of panels). The GRB distributions in the log(l + z) — loge plane are also shown 
in the right group of panels. The symbol style is the same as Figure 1. The contours are the same as that in Figure 1 for all GRBs in our 
sample. 




Fig. 4. — Comparisons of (a) the XRT lightcurves and (b) the Lx i2h ~ ^iso distribution among the low-£ GRBs with Tgo iZ < 5 s (red), 
high-e GRBs with Tgo, z < 2 s (blue), and 44 long-duration high-e GRBs from the Liang et al. (2009) sample (black). Here Lx i2h is the 
X-ray luminosity in the XRT band at the rest-frame 12 hours post the GRB trigger. 



